Negative differential thermal resistance and thermal transistor 
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We report on the first model ol a thermal transistor to control heat flow. Like its electronic 
counterpart, our thermal transistor is a three-terminal device with the important feature that the 
current through the two terminals can be controlled by small changes in the temperature or in 
the current through the third terminal. This control feature allows us to switch the device between 
"off" (insulating) and "on" (conducting) states or to amplify a small current. The thermal transistor 
model is possible because of the negative differential thermal resistance. 

PACS numbers: 44.10.-K, 05.45.-a, 66.70.+f 



The invention of the transistor |l| and other relevant 
devices that control the electric charge flow has led to 
an impressive technological development. In recent years 
some interesting progress has been made also in the con- 
trol of the heat current. Indeed a theoretical model 
of a thermal rectifier has been recently proposedQ in 
which the heat can flow preferentially in one direction. 
More recently, new models of thermal diode have been 
devised[3, which allows to improve the efficiency by 
more than three orders of magnitude. 

In this Letter we propose a model for a thermal transis- 
tor which can control the heat flow in analogy to the usual 
electric transistor for the control of the electric current. 
The reason why such a device can, in principle, operate 
is grounded in the phenomenon of "negative differential 
thermal resistance^ (NDTR)which, as shown below, can 
take place in nonlinear lattices. 

Our model of thermal transistor consists of three seg- 
ments, D, S and G. (see Figure 1(a) for its configura- 
tion). The names D, S and G follow from the ones used 
in a MOSFET (Metal-Oxide-Semiconductor Field-Effect- 
Transistor). Each segment is a Frenkel-Kontorova (FK) 
lattice 0, |i| . Segments S and D are coupled (usually 
weakly) via a spring of constant ki n t, while Segment G, 
the control segment, is coupled to the junction particle O 
between segments S and D via a spring of constant ki n t G ■ 
Temperature Tq is used to control temperature To (at 
the junction O) which determines the heat current from 



D to S. The total Hamiltonian of the model writes: 
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Here Xd,%, %s,i, and xo,i are the particles' displace- 
ments from their equilibrium positions in segment D, 
S and G, respectively. Fixed boundaries conditions are 
taken, i.e., x Dfi = x s .q = x G ,o = 0. 

In the following we will present first the results of our 
numerical simulations, then the theoretical explanation. 
In our numerical simulations the lattices are coupled, at 
their first particles, with heat baths at different temper- 
atures To, Ts and Tq, respectively. We use Langevin 
heat baths and we have checked that our results do not 
depend on the particular heat bath realization (e.g. Nose- 
Hoover heat baths). The local heat flux is defined by 
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FIG. 1: (a) Configuration of the thermal transistor, (b) 
Heat current versus the control temperature Tq- Parameters 
are: T D =0.2, V D =1.0, kn=0.2, fc mt =0.05; T s =0.04, Vs=5, 
ks=0.2, Vg=5, kc~l, kintQ—1. Notice that in a wide re- 
gion both Js and Jd increase when the temperature Tq is 
increased. 



J„ = k(x n (x n — and the local temperature is de- 

fined as T n = The simulations are performed long 

enough to allow the system to reach a non-equilibrium 
steady state where the local heat flux is constant in each 
segment. 

As illustrated in Figures 1-2, the model described by 
Hamiltonian can exhibit different useful functions de- 
pending on the parameters values. 

Thermal Switch: We first demonstrate the "switch" 
function of our transistor, namely we show that the sys- 
tem can act like a good heat conductor or an insulator 
depending on the control temperature Tq. This is shown 
in Figure 1(b), where we plot Jq, Js, and Jd versus Tq. 
When Tq increases from 0.03 to 0.135, both Jd and Js 
increase. In particular, at three points: Tqk 0.04, 0.09 
and 0.135, we have that Jd=Js and thus Jq is exactly 
zero. These three points correspond to "off" , "semi-on" 
and "on" states, at which Jd is 2.4 x 10~ 6 , 1.2 x 10~ 4 and 
2.3 x 10~ 4 , respectively. The ratio of the heat currents at 
the "on" and "off" states is about 100, hence our model 
displays one important function - switch - just like the 
function of a MOSFET used in a digital circuit. 

Thermal modulator/amplifier: As demonstrated 
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FIG. 2: Heat current versus the control temperature Tq. 
Here: T D =0.2, Vd=1.0, k D =0.2, T s =0.04, Vs=5, fc s =0.2, 
femt=0.05, Vg~5, ka=l, fejnt G =0.1. The shadow region is 
the range of variation of Jg in the temperature interval Tg£ 
(0.05,0.135). 



above, the heat current from D to S can be switched 
between different values. However, in many cases, like in 
an analog circuit, one needs to continuously adjust the 
current in a wide range. In Figure 2 we demonstrate this 
"modulator/amplifier" function of our transistor. Here it 
is seen that in the temperature interval TqG (0.05, 0.135), 
the heat current through the segment G remains very 
small (— 10~ 5 ~ 10~ 5 ), within the shadow strip in Figure 
2, while the heat currents Js and Jd are continuously 
controlled from 5 x 10 -5 to 2 x 10~ 4 . 

Let us now turn to the discussion of the underlying 
physical mechanism which allows the transistors func- 
tion. In all cases considered here, the heat resistance of 
segment G has been taken small enough so that Tq ~Tq. 
For fixed thermal baths temperatures Ts and Td, we 
would like now to understand how the heat currents 
Js and Jd depend on temperature To- This depen- 
dence is determined by the differential thermal resistance, 

jd / dJs j and R ( ) Tn 

5 ~ V aTo /T s =const D ~ \ dTo J t d =const 

typical situations the heat current increases with the in- 
crease of the thermal gradient and therefore both Rs and 
Rd are positive. As a consequence, the "current ampli- 
fication factor 1 ^ , 
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namely, in order to induce a change AJd, the control 
heat bath has to provide a larger AJq- This means that 
the "transistor" can never work! 

In order to get the amplification effect, namely a > 1, 
either Rs or Rd must be negative. Thus we must have 
a "negative differential thermal resistance" (NDTR). In 
such case the possibility arises that when To changes 
both Js and Jd change simultaneously in the same way. 
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FIG. 3: (a) Configuration of the model exhibiting the negative 
differential thermal resistance (see definition in the text), (b) 
Heat current versus temperature Tl (at fixed Tr=0.2) for dif- 
ferent coupling constants, kint, with lattice size Nl~Nr=25. 
The system parameters are: Vl—5, Vr—1, kt,=l, kn=0.2. 
(c) Same as (b) but for different system sizes and &i nt =0.05. 
Notice that when Tl<0.1 the heat current increases with de- 
creasing the external temperature difference. 



Therefore the condition Js—Jd can be achieved for sev- 
eral different values of To (as shown in Figure 1). Anal- 
ogously, the condition J$ « Jd can be verified in a wide 
range of To (as shown in Figure 2). In these situations 
heat switch and heat modulator/amplifier are possible. 
In the ideal, limiting case of Rs=—Rd which, in princi- 
ple, can be obtained by adjusting parameters, the tran- 
sistor works perfectly. 

Figure 3 illustrates the mechanism of NDTR. In Fig- 
ure 3(b) we plot the heat current versus temperature Tl 
(see Figure 3(a))for fixed Tr, fixed lattice size and for 
different coupling constants ki nt . In Figure 3(c) instead 
the heat current is plotted for different lattice sizes with 
fixed coupling constants ki n t- 

As it is clearly seen, in a wide range of system pa- 
rameters, for a fixed temperature Tr, there exists a tem- 
perature interval in which a larger value of 7l(< Tr), 
corresponding to a smaller thermal gradient, can induce 



a larger heat current. 

The phenomenon of NDTR can be understood from 
the mismatch between the phonon bands of the two in- 
terface particles. This is illustrated in Figure 4 where we 
plot the power spectra of the two interface particles at dif- 
ferent temperatures Tl for a fixed Tr=0.2 (Figure 3(a)). 
As it is seen, the overlap of the power spectra increases 
as Tl increases i.e. as the thermal gradient decreases. 
Therefore, by increasing Tl we are in the presence of two 
competing effects: from one hand the thermal gradient 
decreases; on the other hand the band overlap increases. 
The beaviour of the heat current depends on which of 
these two effects prevails. In particular the phenomenon 
of NDTR can take place. 

The physical mechanism of such behaviour is due to the 
fact that the variation of the bands position with tem- 
perature depends on the strength of nonlinearity. This 
generalproperty is at the base of the thermal diode mod- 
els 0,0. In particular, for the model discussed in this 
paper, as the temperature Tl increases, the kinetic en- 
ergy of the left interface particle becomes large enough to 
overcome the on site potential barrier and low frequency 
modes appear. 

Finally, we would like to discuss the possibility of build- 
ing a thermal transistor in real physical systems and 
the possible applications. Due to the fast developing 
semiconductor technology, it is now possible to fabri- 
cate several low dimensional devices, such as nanotubes, 
nanowires and thin films in nanometer scale 0- A nan- 
otube with small diameter or a nanowire is basically a 
one dimensional heat conductor. For the purpose of our 
model, which is indeed one dimensional, one can fabricate 
a T-shape nanotube or nanowire. To build the on-side 
potential in the FK model, one may put different seg- 
ments of nanotubes/nanowires on the surface of different 
substrates. 

Among the possible applications we would like to men- 
tion, as an example, that in our model, by adjusting the 
temperature of the third terminal, one can amplify the 
current through the other two terminals. This can be 
used as an efficient heat pump to dissipate heat from a 
microelectronic or nanoelectronic device where a large 
amount of redundant heat is produced due to the fast 
operational speed of the device. This could be very im- 
portant since, as it is known, the redundant heat, if it is 
not dissipated immediately, may deteriorate the device. 

In conclusion, we have built a theoretical model for a 
thermal transistor. The model displays two basic func- 
tions of a transistor: switch and modulator/amplifier. 
The crucial element of the thermal transistor is the neg- 
ative differential thermal resistance. Although at present 
this is no more than an abstract model, we believe that, 
sooner or later, it can be realized in a nanoscale system 
experiment. After all the Frenkel-Kontorova model used 
in our simulations is a very popular model in condensed 
matter physics [5l |6(. For instance, it has been used to 
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FIG. 4: The power spectra of the particles at the right (R) 
side and the left (L) side of the interface, for different tem- 
perature Tl with fixed Tr=0.2. The system parameters are: 
V L =5, Vr=1, k L =l, k R =0.2, k int =0.05, 7V=50. Notice that 
the overlap of the spectra increases as the temperature Tl in- 
creases, which results in the increase of the heat current with 
decreasing the external temperature difference when Tl<0.11. 



model crystal dislocations, epitaxial monolayers on the 
crystal surface, ionic conductors and glassy materials, an 
electron in a quasi-lD metal below the Peierls transi- 
tion, charge density waves, Josephson junctions chains, 
and dry friction. More recently, this model has been em- 
ployed to study transport properties of vortices in easy 
flow channels [8| and strain-mediated interaction of va- 
cancy lines in a pseudomorphic adsorbate system^. 
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